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Abstract 

I calculate the anomalous dimension governing the Q 2 evolution of the gluon (and structure 
functions) coming from the running coupling BFKL equation. This may be expressed in an exact 
analytic form, up to a small ultraviolet renormalon contribution, and hence the corresponding 
splitting function may be determined precisely. Rather surprisingly it is most efficient to expand 
the gluon distribution in powers of a s (Q 2 ) rather than use the traditional expansion where all 
orders of a s ln(l/x) are kept on an equal footing. The anomalous dimension is very different from 
that obtained from the fixed coupling equation, and leads to a powerlike behaviour for the splitting 
function as x — > which is far weaker, i.e. ~ x -0 ' 2 . The NLO corrections to the anomalous 
dimension are rather small, unlike the fixed coupling case, and a stable perturbative expansion is 
obtained. 
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1. Introduction. 

Small x physics has been an active area of theoretical research in the past few years, largely 
due to the first data for x < 0.005 being obtained by the HERA experiments [1] [2]. The crux 
of the debate has been whether the standard DGLAP approach based on renormalization group 
equations and conventionally ordered simply in powers of a s (Q 2 ), or the BFKL equation [3], which 
sums leading logarithms in (l/x), is most effective, and/or whether the two approaches need to 
be combined in some way. While the DGLAP approach has been relatively successful, albeit with 
some significant problems (a valence-like input gluon, undershooting for x ~ 0.01 at the highest Q 2 , 
see [4]), the original BFKL prediction of a behaviour of the form x~ x at small x, with A ~ 0.5, was 
clearly ruled out. A combination of the two approaches, using the BFKL equation to supplement 
the Altarelli-Parisi splitting functions with higher terms of the form a™ +1 ln n (l/x) had some success 
(so long as one avoided factorization scheme ambiguities by working in physical quantities) [5] , but 
this was difficult to sustain with the most recent data. Moreover, the subject was thrown into some 
confusion by the calculation of the NLO correction to the BFKL equation [6] [7] . 

In order to illustrate this I begin with a brief discussion of the LO fixed coupling BFKL 
equation. Working in moment space, i.e. defining the moment of the structure function by 

F(N, Q 2 ) = f x N - 1 F(x,Q 2 )dx, (1.1) 
Jo 

and similarly for the parton distributions (scaled by x), the BFKL equation is 

f(k 2 ,a s /N) = f I (k 2 ,Q 2 ) + ^ T ^K ( q 2 ,k 2 )f(q 2 ), (1.2) 

^ v Jo Q 

where f(k 2 ,a s /N) is the unintegrated gluon four-point function, fi(k 2 ,Ql) is the zeroth order 
input, a s = (3/7r)a s , and the LO kernel is defined by 

Ko(q 2 , k 2 )f{q 2 ) = k 2 ( n f h l~ fi 2 k P + /(fc2) A ■ (1-3) 
V I « -9 I (4(? 4 + fe 4 )2/ 

It is convenient to define the input by fi(k 2 ,Q^) = 5{k 2 — Qq), where in the case of deep-inelastic 
scattering, where one end of the gluon ladder is at a hard scale Q 2 , while the other end is formally 
on-shell, Ql is just a collinear regularization which we let — > ultimately. The "gluon structure 
function" is then given by 

f Q2 rik 2 

g(Q 2 ,N) = j —f(N,k 2 ,Q 2 )xg B {N,Q 2 l (1.4) 

where gB(N, Q$) is a bare, nonperturbative gluon density in the proton which implicitly absorbs 
the collinear divergences in f(k 2 ). The BFKL equation is most easily solved by taking the Mellin 
transformation to 7-space, i.e. 

<»oo 

/( 7 , N) = / dk 2 (k 2 )~ 1 ~" f f(k 2 , N), (1.5) 
Jo 
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where it reduces to 

/(7, N) = fi(7,Ql) + (a s /N) X o(7)f(7,N), (1.6) 
where /(7,Qo) = ex P( — 7^ n (Qo)) an d x(7) * s the characteristic function 

Xo(7) = 2^(1) -#y)~ ^(1-7)- (1.7) 
A little manipulation leads to the expression 

s ^ = h£«^™» wt&^Y (L8> 

This inverse transformation is dominated by the leading pole at 1 — (a s /N)xo(l) = 0, and the 
solution is 

(L9> 

The anomalous dimension j (a s /N) may be transformed to x-space as a power series in a s ln(l/x), 
and has a branch point at N = 41n2a s (at which 7^5) leading to asymptotic small x behaviour 
for the splitting function 

xP°(x) -> 0.07a s x- x /(a s ln(l/x)) 3/2 . (1.10) 

One can, if one ignores the running of the coupling, proceed through exactly the same sort of 
arguments including the NLO correction to the kernel. In this case the "intercept" for the splitting 
function is shifted from A = 41n2a s to A = 41n2a s (l — 6.5a s ). This is clearly a huge correction, 
and implies the breakdown of the perturbative expansion for this quantity. Moreover, the power 
series for the splitting function is dominated by the NLO corrections at all values of x below about 
x = 0.01. For example, using the formulae in [6] the first few terms in the power series for P(x) go 
like 

xP{x) =a s + 2.4q^ 3 /6 + 2.1a^ 5 / 12 + ■ • ■ 

- a s (0.43a s + 1.6a^ + 11.7a 3 £ 2 / 2 + 13.3a^ 3 /6 + 39.7a^ 4 / 24 + 169.4a^ 5 / 12 + •••), 

(1.11) 

where £ = ln(l/x). Clearly, the size of the coefficients more than compensates for the extra power 
of a s . More careful calculations (discussed a little more below), including the running coupling but 
using saddle-point approximations imply the same result, but with the coupling being evaluated 
at Q 2 . This does nothing to aid the convergence, particularly at values of Q 2 ~ lGeV 2 where the 
perturbative analysis of structure function evolution does take place. 

Hence, this NLO correction left the whole question of how to address the evolution of structure 
functions at small x in question, and various alternatives have been proposed to help stabilize 
the calculation. In this paper I demonstrate that the correct way in which to calculate the Q 2 
evolution at small x is indeed to supplement the conventional AP splitting functions with higher 
order corrections obtained from the solution to the BFKL equation, but to take the running of the 
coupling constant fully into account in this equation. The details of this are presented below. 



2 



2. BFKL Equation for Running Coupling. 

Beyond leading order it is impossible to ignore the running of the coupling, since at NLO 
ultraviolet regularization is required, resulting in a ln(/c 2 / fi 2 ^) term where is the renormalization 
scale. Such a term may be eliminated by using the running coupling constant evaluated at the scale 
k 2 . Since this is unavoidably forced upon us at NLO, it seems sensible to consider the fixed coupling 
LO BFKL equation as just a model which would apply in a conformally invariant world, and more 
realistically to work with the BFKL equation with running coupling [8] [9] from the beginning. 
Doing this we obtain 

f(k 2 , Ql a s (k 2 )/N) = fl (k 2 , Ql) + ^ffl ^K (q 2 , k 2 )f(q 2 ), (2.1) 

JV Jo Q 

where 

a s = l/(p ln(k 2 /A 2 )), (2.2) 

[3q = (11 — 2iV//3)/(47r), and Nf is the number of active flavours. One can solve this equation in 
the same type of way as for the fixed coupling case, i.e. take the Mellin transformation with respect 
to (k 2 /A 2 ). It is most convenient first to multiply through by ln(fc 2 /A 2 ), and then obtain 

where /?o = (irf3o/3). Hence, the inclusion of the running coupling has completely changed the form 
of our double Mellin space equation, turning it into a first order differential equation. This has a 
profound effect on the form of the solutions. The equation may easily be solved giving, 

/( 7 ,iV) = exp(-X (7, iV)Z(A)iV))^ 00 dfl ^ Q & exp(X o (7)/0W)<*7, (2.4) 

where 

*o(t) = /^o(7)d7= (2V'(l)(7-|)-ln( T ^^y))- (2.5) 

The leading singularity in the 7 plane for exp(— Xq^/^qN)), is cancelled by an integral from 
— > 7 of the integrand depending on 7 [11], and so up to 0(A 2 /k 2 ) corrections (2.4) simplifies to 

)2\ 



f( 7 ,N) = exp(-X (7)/(A)AO) H dfl{l ' N ' Qo) eMX (l)/(PoN))d^, (2.6) 

Jo "7 



and hence 



1 , 

G(Q 2 ,N) = — I -exp( 7 ln(Q 2 /A 2 ) - X ( 7 )/(A>JV))d7 

J — — ioo T 

x / exp(-7ln(Q 2 /A 2 )+X (7)/(^iV))d7<7B(Q^iV) 
Jo 



2 — IOO 



(2.7) 



= Q E (Q 2 ,N)g I (Q 2 ,N)g B (Q 2 ,N). 
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The essential expression, exp(Xo( / y)/(^oN)), contains singularities at all positive integers, and 
Qi(Qo,N) is not properly denned, since the integrand has singularities lying along the line of 
integration. These are due to the divergence of the coupling at low k 2 and can only be removed by 
some infrared regularization. However, since this factor is independent of Q 2 , it does not contribute 
at all to the evolution of the structure function, and is irrelevant for the evolution. The function 
Qe(Q 2 ,N) is determined by the singularities of exp(— X (^)/(/3oN)) in the 7 plane. This leads to 
a fundamental difference between the cases of the fixed and running couplings. Whereas previously 
the leading singularity was a pole at {a s /N)x{~f) = 1, he. at 7 — > \ as N — ► 41n2a s , now the 
leading singularity is an cut at 7 = 0: there is no powerlike behaviour in Q 2 . Similarly, the branch 
point in the N plane at 41n2a s has become an essential singularity at N = 0: there is no powerlike 
behaviour in x. The introduction of the running of the coupling has changed the character of the 
solution completely. 

Acknowledging that the only real information contained in Ge(N,Q 2 ) is on the evolution of 
the structure function, i.e. defining 

Ge{N,Q ) gives us an entirely perturbative effective anomalous dimension governing the evolution 
of the gluon structure function. The time-honoured technique for solving for Qe{X, Q 2 ) is to expand 
the integrand in (2.7), about the saddlepoint. This results in a contour of integration parallel to 
the imaginary axis, with real part — ► | for the small x solutions. Using this results in an anomalous 
dimension 

00 

T(N, Q 2 ) = l0 {a s {Q 2 )/N) + ^(-/3 a s (Q 2 ))"7n(« s (Q 2 )/iV), (2.9) 

n=l 

i.e., the effective anomalous dimension is the naive leading order result with coupling at scale Q 2 
plus a series of corrections in increasing powers of — (3oa s (Q 2 ). However, each of the ^(a s (Q 2 )/N) 
is singular at N = X(Q 2 ), and the power of the singularity increases with increasing n [12]. Hence, 
although the series for the resulting splitting function is in the small quantity a s (Q 2 )(3o, the ac- 
companying coefficients are progressively more singular as x — > 0. The saddlepoint approximation 
is therefore not a reliable result as x — > and explicit investigation reveals that it is only really 
quantitatively useful when a s (Q 2 ) ln(l/x) is so small that the effective anomalous dimension is 
effectively the LO in a s part xP(x) = a s (Q 2 ). Therefore any calculations of the anomalous di- 
mension which rely on an expansion about the saddle-point lead to very inaccurate and misleading 
results for small x. 

This instability is not surprising if one examines the integrand along the saddle-point contour 
of integration, noting that it is very different from the Gaussian form the saddle-point method 
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assumes [12], and also if one notes that it is an expansion obtained from approaching 7 = | and in 
terms of functions of N which are singular at N = X(Q 2 ), whereas we know that the full solution 
no longer sees these points as anything special. In fact, the knowledge of the singularity structure 
of the integrand implies that 7 = is a more fruitful point on which to concentrate. Prompted by 
this we may move the contour of integration to the left and simultaneously use the property that 
the integrand dies away very quickly at infinity to close the contour so that it simply encloses the 
real axis for 7 < 0. It is then useful to express Xo{l) m the form 

00 

X0 (7) = 1/7 + E 2 ^ 2n + ( 2 - 10 ) 

71=1 

which is strictly valid only for | — y | < 1. Doing this we may write 

*o(7) = Hi) +7E + j>2 2^ T t T7 2 " +1 , (2.11) 

71 = 1 

and the integrand for Ge(N,Q 2 ) becomes 

7 -i/(**)-i exp L t _ -J—( 1E + f; a n7 2n+1 )) , (2.12) 

where t = ln(Q 2 /A 2 ) and a n = 2((2n + l)/(2n + 1). The contribution to the integral from 
— > — 00 + it is now the same as that from — 00 — ie — > up to a phase factor, and we may write 

(2.13) 

where the integral has to be understood as an analytic continuation, since there are singularities 
along the real axis, and strictly speaking the integrand is well defined only for 7 > — 1. Since the 
factor of exp(7i), is present this latter point leads to an ambiguity of order exp(— t), i.e. C(A 2 /Q 2 ) 
into the value of GE{N,t). 

In order to evaluate the above integral it is convenient to let y = yt, resulting in 

,0 / , 00 \ (2.14) 

The latter exponential may be expanded as a power series in y/t and the integral evaluated using 
the standard result that 

(_l)T(-l/(A)iV) +n)= I y-Vtfo*)-! exp(y)y n d 7 , (2.15) 

J —00 
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and hence 

g E (N, t) = - sm(^y) exp ("T^fr) ) t 1/W ° N) £ A n (l/ (fa N))t~ n (-1)^-1/ (PoN) + n). 

(2.16) 

This result was noted in [13], as was the fact that it may be simplified by using the relationship 
that as N -> 0, F(-l/(p N)+n) -> r(-l/((3 N))(-l/(i3 N)) n . However, it is important to notice 
the more general result that for all N 

(_1)T(-1/(A,A0 + n) = r(-l/(A,JV))A n (-l/(A,iV)), (2.17) 

where 

n 

A n (-1/(P N)) = (-1)" J2 (-l) m dnm(PoN)- m , (2.18) 

m=l 

where d nm are positive coefficients and d nn = 1. Hence, ignoring the common factor of 
— sin(7r/(/3 iV))r(— l/(p N)) exp(— 7e/( / 5 A r )), which has no i dependence, and is irrelevant for 
the anomalous dimension, 

G E (N,t) = t 1 /^ (l + f2A n (l/(p N))t- n A n (-l/(P N))) (2.19) 

^ n=3 ' 

where the A n are simply calculable from the expansion of ex.p^—l/(P( ) N)^2'^' =1 a n (y/t) 2n+1 ^j. 

The common factor of t 1 /^ ^) is the well-known double-leading-log result coming from just the 
LO a s (Q 2 )/N part of the anomalous dimension. Multiplying this we have an expansion as a power 
series in 1/t or equivalently in a s (Q 2 ). In fact 

n 

t- n A n (-l/(p N)) = (a s (Q 2 )/N) n d nm {-(5 a s {Q 2 )) n - m {a s {Q 2 )/N) m - n . (2.20) 

m=l 

This explicitly demonstrates that we obtain a set of running coupling corrections to a LO result. 
Substituting this expression for Ge{N,£) in (2.8) one obtains an expression for the anomalous 
dimension as a power series in a s (Q 2 ), where at each order we have the leading divergence in 1/N 
plus a sum of running coupling correction type terms. With a little work one may regain the whole 
leading j (a s (Q 2 )/N) (though it is necessary to keep some subleading terms in the A n to do this), 
along with a tower of terms which are subleading in powers of (3oa a {Q 2 ) to this leading anomalous 
dimension, i.e. one obtains all the corrections to this naive LO anomalous dimension due to the 
running of the coupling. 

The general features of this full, running coupling BFKL anomalous dimension may be appre- 
ciated quite easily. The important fact to note is that although the A n (—1/((5 N))) — > (l/{(3oN)) n 
as N — > 0, the function oscillates with l/(/3 A r ), and remains very much smaller in magnitude than 
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this asymptotic form until very small N, roughly until l/N > n. This coupled with the accom- 
panying factor of t~ n means that for reasonable t, i.e t > 5 (Q 2 ^ lGeV 2 ), only the first 5 or so 
terms in (2.19) contribute for N > 0.25. Hence, to a very good approximation 

g E (N,t) _,./<*».(, _ _||_ A 3 ( _ 1/( a W) ) _ _||_ A5( _ 1/(/5oW)) ), (2 . 21) 

and in fact the smallness of the coefficient makes even the t~ 5 term almost negligible. GE(N,t) 
initially grows as N falls due to the t 1 /(^° Ar ) term. However, for N ~ 0.6 the negative contribution 
from the t~ 3 term starts to contribute and ultimately drives the gluon structure function to negative 
values. The result is shown in fig. 1. dQ e(N ,t) j dt may simply be evaluated also using (2.21), and 
shows the same general shape, but does not become negative until a slightly lower value of iV as 
also seen in fig. 1. Hence the anomalous dimension develops a pole at a finite value of N, given 
quite accurately by 

JC(3W^ + (2 22) 

where we use the explicit form of A 3 (— l/(/3oA r )). The value of N for this leading pole is shown 
as a function of t in fig. 2, and for the sort of values relevant at HERA is ~ 0.25. Going to 
N < 0.25 higher order terms become important, and the positive l/((f3oN) 2 t 6 )A 6 (— \/([3qN)) term 
pulls C/E(N,t) back to positive values, and another pole, with opposite sign residue, appears in 
T(N,t). At even lower N the analytic expression eventually breaks down, but numerical results 
show a series of poles becoming closer together. Nevertheless, the position of the leading pole is 
essentially determined by the first handful of terms in the power series in a s (Q 2 ) for QE(N,t), and 
hence so is the asymptotic behaviour of the small x splitting function, i.e. ~ x -0 ' 25 . Hence, the 
introduction of the running coupling has a dramatic effect on the singularity structure of the LO 
BFKL anomalous dimension, turning the cut into a series of poles, and changing the position of 
the rightmost singularity by a factor of ~ 0.4. This result of the pole in the anomalous dimension 
has been already noted in [14] using numerical techniques, and in the context of a resummed NLO 
anomalous dimension. Here I particularly stress the huge modification of the naive LO BFKL 
anomalous dimension. This huge change is apparent over a wide range of N, and in fig. 3. a I show 
the anomalous dimension as a function of N for all values right of the leading singularity. As one 
sees, it is rather closer to the simple a s (Q 2 )/N expression than to the naive BFKL result. 

Here I should comment on the limit of the analytic expression. As noted, it involves a series 
expansion not valid over the whole contour of integration. This is reflected in the fact that the 
overall magnitude of the A n (—l/(f3oN)) increases like n! in general. This means that the series in 
(2.21) is actually asymptotic, although it is an oscillating series, so in principle is unambiguously 
resummable (i.e. the integral in (2.7) does exist). However, the greatest accuracy may be obtained 
from (2.21) by truncating the series at order n ~ t. In practice I always use n = 5. Substituting 
the truncated expression for 0E(N,t) then results in an infinite series in a s (Q 2 ) for T(N,t) which 
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is convergent for any given N right of the leading pole. The accuracy of the analytic expression 
can be found by comparing with results obtained from evaluating (2.7) using numerical integration, 
and for the gluon structure function for A to the right of the leading pole is found to be rather 
better than 0.1% for t = 6 and falls like exp(— t). 

In order to investigate the quantitative effect of the BFKL anomalous dimension on structure 
function evolution it is necessary to calculate the BFKL splitting function as a function of x. This 
is where an analytic expression for the anomalous dimension is particularly useful. A series of 
numerically obtained values of T(N,t) allows an approximate determination of P(x,t), but it is 
very difficult to be accurate, especially for the wildly oscillating functions of 1/N which do in fact 
make up Ge(N, t). However, I now have an explicit series for T(N, t) in powers of a s (Q 2 ), obtained 
from the truncated expression for Qe{ji^€). The A-dependent functions at each power of a s (Q 2 ) 
of course become larger at small A as the series progresses, and to reach small enough x more and 
more terms are needed. However, at a fixed value of A there is no such growth, and the same 
is therefore true for fixed x. Hence, one only needs to work to a finite order. Limiting oneself 
to x > 10~ 5 and t > 4.5, the suppression of the A n (— l/(/3 A)) is quite significant and seventh 
order in a s {Q 2 ) is easily sufficient. The splitting function for t = 6 is shown in fig. 4. a. One sees 
that it is hugely suppressed compared with the naive LO BFKL splitting function, and is even 
lower than the 0{a s {Q 2 )) contribution for x between about 0.1 and 0.001. The fact that there 
is deviation from the standard NLO in a s (Q 2 ) splitting functions tells us that BFKL influenced 
structure functions are important. The ultraviolet renormalon contribution is approximated by 
constructing an x-space function which matches the iV-space results for a variety of TV. Although 
this contribution turns out to be a larger fraction of the total than in A-space, it still only makes 
a very small correction to the evolution. This will be discussed more in a subsequent paper [15]. 

3. Small x Structure Functions. 

It is important to realize that all of the above results are in a sense ambiguous because they 
deal with a particular way of defining the gluon parton distribution. It is defined in the natural way 
for a discussion of BFKL physics, but nonetheless is really an intrinsically factorization-scheme- 
dependent quantity. One may define a real structure function by including a hard scattering cross 
section at the top of the gluon ladder. This modifies (1.4) to 

F l (Q 2 ,N) = a s ^ -^a^k 2 /Q 2 )f(N,k 2 ,Q 2 )g B (N,Q 2 ). (3.1) 

Pulling out the overall factor of a s and taking the Mellin transformation of (3.1) with respect to 
(Q 2 /A 2 ) leads to the simple expression 



F i { 1 ,N) = h i , g { 1 )g{ 1 ,N). 
8 



(3.2) 



Thus we may solve for Fi{N, t) in the same way as for Q(N,t) obtaining the same divergent Q 2 - 
independent part and a Q 2 -dependent part given by solving 

T ^(N,t) = —J i -^nl e Mlt-X (l)/WoN))dj. (3.3) 

We may proceed as with the gluon structure function by expanding the ^,(,(7) (which were calcu- 
lated in [16]) as a power series about 7 = 0. This results in an expression 

F E ,(N, t) = i^oAO fi + jr B hn (l/(P N))t- n A n (-l/(P N))) , (3.4) 

^ n=l ' 

where the B itn (l/((3oN)) are now determined not only by the power series in 7 obtained from the 
expansion of A (7), but also from the expansion of h itg (^). In particular they now contain parts 
at zeroth order in l/(p N). In order to draw the most direct analogy to the gluon we define 

All the terms in the expression for this physical anomalous dimension, determined from (3.4), are 
part of the standard LO T® LL (a s {Q 2 ) / N) or are subleading by powers of Poa s (Q 2 ) to this, i.e. again 
we obtain the naive LO result plus running coupling induced corrections to this, and the "coefficient 
function" ^,3(7) contributes only to the running coupling corrections. 

TLL(N,t) is shown along with the gluon anomalous dimension in fig. 3. a. Clearly the effect of 
the additional coefficient function, and hence additional running coupling corrections, is to make 
^LL(N,t) dip significantly below the simple a s (Q 2 )/N and to reduce the value of the intercept 
compared to the gluon structure function. This is reflected in the effective splitting function 
PhL{x,t) which is shown in fig. 4. a. This time the dip below the 0(a s (Q 2 )) part is far more 
pronounced. Also, going to low enough x, we see that the splitting function turns over again, 
showing that the subleading poles in the anomalous dimension may have large residues compared 
to the leading pole, and the increase in Pll{x) with decreasing x is not monotonic. This illustrates 
that as far as phenomenology at HERA, or any forseeable collider, is concerned, the value of the 
intercept for the anomalous dimension is simply not relevant to the evolution of structure functions. 
I should note that in the case of the physical anomalous dimension and splitting function the power 
corrections due to the nonconvergence of the series are somewhat larger, about 1% in A-space for 
t = 6, and lead to a small, but not insignificant effect on the evolution in x-space, and must be 
accounted for. This will be discussed more in [15]. 

One can follow exactly the same procedure for the other physical anomalous dimension 

9 ^ Q f ] =r 2 UQ 2 ,N)^ L (N,Q 2 ), (3.6) 

obtaining qualitatively very similar results. Again the splitting function initially dips as x decreases, 
but grows again in the same way as Pll(x, t) at lower x. This has extremely important implications 
for the evolution of F 2 (x, Q 2 ) in the HERA range, and this will be presented in detail in [15]. 
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4. NLO Corrections. 

So far I have demonstrated that using a s (k 2 ) in the BFKL equation, as in (2.1), has a profound 
effect on the form of the solution for the anomalous dimension. However, it is necessary to check 
that the results presented are not severely modified by the inclusion of the NLO kernel, i.e. the 
perturbative calculations are stable, and also to justify that the choice of scale in the coupling is 
correct, or at least leads to accurate results. The NLO kernel was presented in [6] and the way in 
which to solve at NLO with a running coupling was presented in [17]. Writing the NLO equation 
as 

f{k 2 , Ql) = fi(k 2 , Qg) + {^P-) [ ^r(Ko(q 2 , k 2 ) - a s (k 2 )K 1 (q 2 , k 2 ))f(q 2 ), (4.1) 

using just the one-loop coupling leads to a 2nd order differential equation in 7-space 

t&bm - fllhM _ J_ d(Xo(7)/(7,*0) _ * l{l) f {l N) (4 2) 

d 7 2 d 7 2 p N d 1 3P 2 N X1[1)J[1, h { ] 

This can be solved in a very similar way to LO, i.e. it factorizes into the same form as (2.7) with 
Q 2 -dependent part given by 

Q E ,NLo(N,t) = -L j 2+ Iexp( 7 t-X JVLO (7,Af)/(^ Af)) a !7. (4.3) 



2 ' 



However, X^Loil-,^) is rather more complicated than the previous -^0(7). It can be expressed in 
the form 

X NLO (7,N)= I xnlo ( A /,N)dj, (4.4) 

where XNLoil-, X) can be written as a power series in iV beginning at zeroth order with Xo(7)- As 
seen in [17], though here ignoring any resummations in N, the explicit form is 

*™(7, ») - X.M - + ^ (- (SM)' - A (*M)' ) + . . . , (4.5) 



Xo(7) Xo V VXo(7)y VXo(7) 

where ^2(7) would also appear at order iV 2 , had I included it. I shall generally ignore all but the 
first two terms. 

Firstly, I shall address the choice of scale. It was known in [18] that the correct scale was really 
(k — q) 2 , but that k 2 could be used, leading to a part of the NLO kernel which is proportional to 
(3q, i.e. there is a contribution to Xi( 7 ) of the form |A)(Xo(7) + Xo(7))- Substituting this into 
(4.3) leads to a contribution in the integrand of the form exp(i(ln(xo(7)) + -^0(7)))- This can be 
expressed as a power series which at low orders is 1 + I.67 3 + I.247 5 . Hence, this scale-induced 
factor has the same form as the h^gi^y), and not surprisingly results in additional running coupling 
corrections to the anomalous dimensions. However, the terms in the series do not start until third 
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order, have small coefficients, and have an effect much smaller than ^,l(2) (t)- Hence, the correction 
for this "incorrect" choice of scale is very small, though in principle it seems as though the factor 
just considered should really be taken as part of the LO result, since it just gives running coupling 
corrections only. 1 

So the choice of a s (k 2 ) is in practice very reliable, and may easily be corrected for. We must 
now consider the rest of the NLO correction to the kernel, which is much larger. Here we have 
an ambiguity in precisely what the NLO calculation means. Do we simply solve (4.1), producing 
the infinite series in (4.5)? Do we truncate Xnlo(i, N) after the second term, and if so do we 

use the whole of exp(l//? J7 (Xi(7)/Xo(7))^7) or j ust expand it out to first order? There are 

2 

particular problems associated with all choices. I choose the NLO definition such that the anomalous 
dimension receives only corrections which are one power of a s {Q 2 ) down on the leading order one, 
i.e. the LO anomalous dimension is of the form T°(a s (Q 2 )/N, p a s (Q 2 )) and the NLO corrected 
one is of the form T°(a s (Q 2 )/N, f3 Q a s {Q 2 )) + a^Q^V^a^Q^/N, f3 a s (Q 2 )). Roughly speaking, 
this involves keeping only the first two terms in (4.5), and expanding exp(l//3 Ji{Xi(.l)/Xo(l))dj), 

2 

out to just first order in l//3 . However, the part of (xi(7)/Xo(7)) behaving like I/7 must be treated 
very carefully. Details will be presented in [15]. 

Solving for the NLO anomalous dimension using the same techniques as at LO results in 
the NLO correction. Unlike the case of fixed coupling, or the naive results of the saddle-point 
evaluation, these corrections are rather small. The positions of the leading poles in the anomalous 
dimensions are shown in fig. 2, and one can see that they change from about 0.21 for Fll at LO 
to 0.16 at NLO, and that the Q 2 -dependence reduces a little. Similarly the anomalous dimension 
t) over a wide range of N shows only a very small change going from LO to NLO. (I use 
the physical anomalous dimensions to avoid any ambiguity. The results are very similar for the 
gluon distribution.) This is shown in fig. 3.b where the part of the NLO anomalous dimension 
at first order in a s (Q 2 ), i.e — 0.935a s (Q 2 ), is not included, since this should properly be included 
at LO in a combined leading order in a s (Q 2 ) and a s (Q 2 ) ln(l/x) expansion scheme. Alternative 
definitions of NLO lead to very similar results except at very high values of N. Note that the NLO 
correction is actually positive for N ~ 0.9 - very different from the case where running coupling 
corrections are not included and all NLO corrections are negative. I should also note that the 
ultraviolet renormalon contribution is a far larger proportion of the NLO correction than it is of 
the LO contribution, being typically 10% for t = 6, and needs to be accounted for [15]. 

One can also make the transformation to x-space and calculate the NLO corrected splitting 
function. This is shown for t = 6 in fig. 4.b, where the contributions oc 5(1 — x) both from the 
0(a s (Q 2 )) part and the running coupling corrections to this are absent. The latter of these is a 

1 In fact, the anomalous dimensions and splitting functions already presented in this paper contain 
corrections from this factor, though this makes only a very small difference. 
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very small contribution. The NLO corrected splitting function is not too different from that at LO, 
as one can see. However, the real import of the NLO corrections as far as physics is concerned is 
the effect it has on the evolution of the structure function. This is demonstrated in fig. 5 where 
the evolution of a suitable model for the structure function Fl(x, Q 2 ), i.e. (1 — x) 6 x~ - 2 , is shown 
both for the LO running coupling splitting function, and for the NLO corrected one (all 5(1 — x) 
contributions other that at first order in a s (Q 2 ) one are included). As one sees, at this (rather 
low) value of t, i.e. Q 2 ~ 8GeV 2 , the effect of the NLO corrections is only of order 10%, and is 
positive for x ~ 0.01. Also shown for comparison is the contribution from the a s (Q 2 )/x term alone. 
The running coupling BFKL splitting function leads to slightly quicker evolution than this latter 
contribution for x *5 0.01, but for x ~ 0.0001 the evolution is significantly suppressed. 

Hence, the NLO corrections to the running coupling BFKL derived splitting function are well 
under control, both in terms of the asymptotic powerlike behaviour of the splitting functions and 
in terms of the evolution in the range currently accessible to experiments. Beyond the running 
coupling corrections no further resummations are necessary, or even useful. This is very much in 
contrast to the case where both ends of the gluon ladder are associated with a hard scale, the 
so-called "single scale" processes. In this shown in [19] and developed in [17] and [14], 

the conventional BFKL expansion is fundamentally flawed due to high order poles at 7 = and 
1, which need to be resummed. Without resummation, all calculations are badly behaved over the 
full range of N, not just at low N. In the case of deep inelastic scattering the collinear factorization 
procedure automatically orders the poles at 7 = correctly, and the above problem shows up in 
high order poles at 7 = 1 only. The anomalous dimension is totally dominated by the region very 
close to 7 = 0, as this paper shows, and is very insensitive to effects at 7 = 1. Including the type of 
resummation in [19] [17] alters results from the NLO corrected case by a very small amount, and is 
possibly no more influential than the remaining NNLO effects for which it does not account. The 
corrected treatment at 7 = 1 is essential if one is attempting to obtain information about the input 
form of the gluon, i.e. Gi{Qq, N), but this, along with the whole subject of single-scale processes, 
is also plagued by the infrared divergence problem. A discussion of such issues can be found in 
[14] [20]. 

5. Conclusions. 

I have shown that an analytic expression for the anomalous dimensions and splitting functions 
obtained from the running coupling BFKL equation may be obtained by a power series solution 
for the Q 2 -dependent part of the gluon structure function in terms of a s (Q 2 ). This has extremely 
good accuracy, with only very small errors, which may be interpreted as ultraviolet renormalon 
contributions, and may be calculated numerically. Moreover, I find the remarkable result that the 
form of this anomalous dimension and splitting function is almost completely determined by only 
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the first handful (~ 5) of terms in the expansion for the gluon. This is in complete contrast with 
the case of fixed coupling, where an all orders summation is needed, and it would be interesting to 
understand the origin of this phenomena. 

My results also prove that the effect of the running of the coupling is to weaken the asymptotic 
powerlike growth of the splitting functions severely compared to the naive results, and even to lower 
the splitting function below the a s (Q 2 )/x contribution for 0.001 ^ x <S 0.2. It also makes the NLO 
correction to the splitting functions relatively small, both for the value of the intercept and for the 
evolution of structure functions for x > 10 -5 , and therefore stabilizes the perturbative series. I also 
note that a previous conjecture that the effect of the running coupling in the BFKL equation could 
be accounted for using an x-dependent scale for the coupling [12], resulting in falling coupling for 
decreasing x, turns out to be largely correct so long as the change in the scale of the coupling is 
moderate compared to the scale itself, though it fails when this criterion is not satisfied. In practice 
this condition is identical to that specifying that diffusion in the fixed coupling BFKL equation 
is not too large and therefore that the virtualities sampled in the running coupling equation are 
not too far away from Q 2 . This results in the requirement that t 3 < 201n(l/x) [21]. This covers 
most of the HERA range for Q 2 > lGeV 2 , and means that a phenomenological analysis using 
the explicit resummation in this paper leads to very similar results to the approach in [12]. In 
particular, the input F^(x,Qq) is of the same shape as F 2 (x,Qq) for as low as lGeV 2 , in 
contrast to the standard NLO-'m-a s (Q 2 ) approach, and the quality of the fit, particularly for small 
x data, is far better than using the conventional approach. A more detailed presentation of such 
a phenomenological study will appear. However, I can say that the inclusion of corrections to the 
fixed order in a s (Q 2 ) splitting functions from the running coupling BFKL equation not only leads 
to a stable perturbative expansion but also to a clear improvement in the comparison to data. 
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Figure Captions 

Fig. 1. The Q 2 -dependent part of the gluon structure function, C/E(N,t), and cLQe{N ,t) / dt 
as a function of N for t = 6 (Q 2 ~ 8GeV 2 ). The (^-independent factor of 
— sin(7r/(/3o-ZV))r(— 1/((3qN)) exp(— 7s/(/3 A r )) is included in both in order to produce 
a more smoother iV-dependent normalization of the functions. 

Fig. 2. The positions of the leading poles in the anomalous dimensions for the gluon structure 
function at LO, for Fl at LO and for Fl at NLO. 

Fig. 3. a. The anomalous dimensions for the gluon structure function at LO and for Fl at 
LO plotted as functions of JV for t = 6 (Q 2 ~ 8GeV 2 ). Also shown is the 0(a s (Q 2 )) 
contribution a s (Q 2 )/N. b. The anomalous dimensions for for Fl at LO and at NLO 
plotted as functions of N for t = 6. 

Fig. 4. a. The splitting functions xP gg (x) and xPll{x) plotted as a function of x for t = 6 
(Q 2 ~ 8GeV 2 ). Also shown is the 0(a s (Q 2 )) contribution a s (Q 2 ), and the naive LO 
BFKL splitting function with coupling a s (Q 2 ). b. The splitting functions xP££(x) and 
xP L y L LO (x) plotted as a function of x for t = 6. 

Fig. 5. The values of dF L (x, Q 2 )/dlnQ 2 , for F L (x,Q 2 ) = x"°- 2 (l - xf, due to the LO splitting 
functions Pll( x ) aR d the NLO splitting function P^ L LO {x) 1 plotted as a function of x 
for t = 6 (Q 2 ~ 8GeV 2 ). Also shown is the evolution due to the 0(a s (Q 2 )) contribution 
P(x) = a s (Q 2 )/x. 
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